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Confexf. Small amplitude oscillations are observed in red giant branch(RGB) stars. Data on such oscillations are a source of infor- 
mation about the objects, notably about properties of convection in their envelopes and about the systems these objects inhabit. The 
OGLE-III catalog contains data for about 80 thousand small amplitude variable red giants (OSARGs) in the Large Magellanic Cloud. 
Aims. We want to explain variability in OSARGs as the solar-like oscillation and to associate the peaks in power spectra with fre- 
quencies of acoustic modes. 

Methods. We use data on reddening-free magnitudes of the objects and interpret them in terms of stellar physical parameters using 
tabulated isochrones calculated for ages and composition parameters corresponding to the upper RGB of the LMC. Massive data on 
the peak frequencies and amplitudes are compared with expectations for stochastically excited oscillations. The frequencies are also 
compared with those calculated for radial modes in envelope models with parameters taken from the isochrones. 
Results. In stars close to the tip of the RGB, the peaks in power spectra are found in the 0.1-1.0 /jHz range, which is consistent with 
extrapolation of the frequency-luminosity relation for the solar-like oscillation. The dominant peaks occur close to the first two radial 
overtones. The increase in amplitude with luminosity is slower than linear. The exponent s ~ 0.9 is similar to what is found from 
recent analysis of CoRoT data on less luminous red giants. Frequency separations between dominant peaks are found to be smaller 
by about 20 % than calculated separations between these modes. After examining various possibilities, we left this discrepancy unex- 
plained. 

Conclusions. The small amplitude variability of stars at the RGB tip is likely to be caused by a stochastic excitation of acoustic 
oscillations but interpreting of individual peaks in power spectra presents a problem. 

Key words, asteroseismology - stars: evolution -stars: late type -stars: oscillations - galaxies: Magellanic Clouds 
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1. Introduction 

It is now generally accepted that acoustic oscillations in the Sun 
are stochastically excited by turbulent convection. Oscillations 
driven in this way, called solar-like oscillations (SLOs), must be 
excited in all stars having an outer convective zone but, owing 
to small amplitudes, they are not easy to observe and study. The 
detection of SLOs in stars along the red giant branch (RGB) was 
among the most interesting results of asteroseismology with the 
space telescopes CoRoT (De Ridder et al. 2009) and KEPLER 
(Stello et al 2010). Oscillations of this type probe stellar convec- 
tion, so it is important to detect them in objects over a possibly 
wide range of surface parameters. Data on mode frequencies in 
individual red giants could give us important constraints on evo- 
lutionary models. This requires, however, a credible mode iden- 
tification. 

The signature of SLOs is an excess of power around a given 
frequency v max , which decreases with growing stellar luminosity. 
The L decrease was expected (Kjeldsen & Bedding 1995) and 
partially confirmed by observations. Data from CoRoT (Mosser 
et al. 2010) and KEPLER (Stello et al. 2010) reach down to 
Vmax ~ 2 //Hz, which roughly corresponds to L = 300L o . This 
number is about one order of magnitude above the red clump 
composed of helium-burning stars and one order of magnitude 
below the RGB tip (TRGB), which stars reach at the moment of 
helium ignition. 

On a large scale, small amplitude oscillations in luminous 
RGB and asymptotic giant branch (AGB) stars were detected as 



a byproduct of microlensing surveys MACHO and OGLE. Using 
data from the first survey, Wood et al. (1999) and Wood (2000) 
show that the LMC red giants form five ridges (sequences) in the 
period-luminosity plane. The ones labeled A and B were asso- 
ciated with the small amplitude pulsators. Later, the MACHO 
data were analyzed by Fraser et al. (2005, 2008). Data from 
the OGLE-II survey for both Magellanic Clouds were the sub- 
ject of independent analyses by Kiss & Bedding (2003, 2004), 
Soszyriski et al. (2004), and Ita et al. (2004a,b). Additional 
period-luminosity sequences formed by low amplitude objects 
were identified by these authors. Interpretations of these se- 
quences in terms of radial pulsation modes in RGB and AGB 
stars were considered. Evidence for the multiperiodic character 
of pulsation was found in both sets of data. 

Here we use the acronym OSARGs OGLE small amplitude 
red giants, which was introduced by Wray et al. (2004) in their 
paper on variability of red giants in the Galactic bulge. Soszyriski 
et al. (2004), who discovered a large number of OSARGs in 
the Magellanic Clouds (15400 in LMC and 3000 in SMC), 
showed that they fall into seven detached ridges in the period- 
luminosity plane. Four of the ridges extending beyond the TRGB 
were attributed to the AGB stars. These objects were named 
type a OSARGs. The ones lying within the three ridges ending 
at TRGB were regarded genuine RGB stars and named type b 
OSARGs. 

The OGLE-III Catalog of Variable Stars (Soszyriski et al. 
2009) contains data for about 8 x 10 4 OSARGs in the Large 
Magellanic Cloud (LMC) and is based on photometric data col- 
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Fig. 1. A sample of power spectra of a and b type OSARGs. For each spectrum we give the OGLE identification number, the 
Wesenheit index (Wi, see Eq. 1), and the amplitude of the highest peak in the the I-band (A/). 



lected during 8 years (2001-2009) of continuous observations 
within the framework of the OGLE-III project. For some of 
the objects, the time base was increased to 13 years thanks to 
photometry from the OGLE-II project. In this paper, we fo- 
cus mainly on type b, because we only have adequate models 
for RGB stars. First, in a phenomenological way, we address 
the question how the OSARGs' variability is driven. Next, we 
present an attempt to explain the dominant peaks in terms of 
stellar oscillation modes. 



2. Are the OSARGs solar-like pulsators? 

The origin of OSARGs is not a new question. It has been already 
addressed by Soszyriski et al. (2007), but some arguments in fa- 
vor of such an interpretation were presented in a relatively short 
section, which went unnoticed by people working on SLOs. 

Stochastic excitation must take place in red giants, but it is 
not at all clear whether the observed variability of OSARGs 
is predominantly caused by this mechanism. In fact, Xiong & 
Deng (2007) find that some low-order radial modes are unsta- 



ble in RGB stars. This instability arises mainly from a small 
excess of the driving effect of turbulent pressure over damp- 
ing effect of the turbulent viscosity. Both effects originate in 
the perturbed Raynolds stress, and are difficult to estimate. Thus 
both the stochastic excitation and the self-excitation of unsta- 
ble modes must be regarded as viable suggestions. In the wake 
of uncertainties in the description of convection in red giants, 
we propose a closer look at observational data to see how the 
properties of low-amplitude oscillations in luminous red giants 
follow the trends seen in fainter objects. 

We began by comparing patterns in power spectra. Figure 1 
shows examples of spectra for OSARGs, based on the OGLE I- 
band photometry. We chose the ones that have prominent peaks 
contributing to a.2, &3, i>2, ^>3 sequences (Soszyriski et al. 2007), 
which are the most prominent in the PL-plane. The patterns seen 
in this figure are reminiscent of those shown in Fig. 1 of De 
Rider et al. (2009) or Fig. 2 of Stello et al. (2010), except that 
the frequency range is a factor 100 lower. In all periodograms 
shown in our Fig. 1, we see the two prominent peaks that were 
the feature behind our selection rule, but besides this, there is 
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Table 1. Selected isochrones and stellar parameters at Wj — 12 - 



Age [Gy] 


[M/H] 


Z 


Y 


M/M 


L/L 




7 


-0.659 


0.004 


0.251 


0.965 


1795 


0.657 


4 


-0.659 


0.004 


0.251 


1.145 


1945 


0.658 


4 


-0.953 


0.008 


0.256 


1.222 


1563 


0.642 



considerable diversity, especially in the amount of power around 
the dominant peaks and in the presence of the secondary peaks. 
In all cases there is a power localized in some intermediate 
frequency range that moves leftwards with increasing luminos- 
ity. Simultaneously, the separation between dominant peaks de- 
creases. These are features seen in all power spectra of SLOs. 

The low-frequency power in the OGLE data may be at- 
tributed to the granulation noise, like in other cool stars. 
However, the high-amplitude peaks in this range come from 
well-known but still unexplained phenomenon of the long sec- 
ondary period (LSP) observed in a significant fraction (~30%) 
of red giants (see, e.g., Nicholls et al. 2009). The power spectra 
shown in this figure are very similar to those shown by Tabur et 
al (2010) for bright Galactic M giants. Some of those objects are 
evidently located near TRGB. 

We turn now to data on all OSARGs that contribute to the 
sequences \>2 and b3. From the OGLE-III Catalog, we took the 
respective frequencies and the Wesenheit index, Wj, for individ- 
ual objects. This index, which is the reddening-free stellar mag- 
nitude, is defined as 



Wj =7-1.55(7-7), 



(1) 



where V and I denote mean magnitudes in the respective bands. 
Points in Fig. 2 present data on the most significant periodici- 
ties (S IN > 5). Up to two points are plotted per star. To com- 
pare these data with information on SLOs in less luminous stars, 
we need stellar parameters, such as luminosity, effective temper- 
ature, and mass. To this end, we used isochrones downloaded 
from the BaSTI Library (Pietrinferni et al. 2006) for ages and 
metal abundance parameters appropriate to stars at TRGB in 
the LMC (Salaris & Girardi 2005). For the distance modulus 
to LMC, we adopted 18.5 mag, which is close to the mean value 
from recent determinations by various methods (Schaefer 2008). 
The characteristics of the selected isochrones and model param- 
eters at the point where Wj = 12 mag are listed in Table 1. In 
the range covered by the data, there is about factor 4 increase in 
luminosity between the highest and the lowest frequencies. 

The three lines in the upper panel of Fig. 2 were calculated 
along the isochrones with the Kjeldsen & Bedding (1995) ex- 
pression 



L Q M ( T \ 3 5 



(2) 



for the frequency corresponding to maximum of the acoustic 
power. This expression is consistent with data on solar-like oscil- 
lation from main sequence up to lower RGB (Bedding & Kjelsen 
2003, Stello et al. 2007). We may see that it also approximately 
applies to our data, which extend up to TRGB. This is an ar- 
gument in favor of the solar-like nature of oscillations in type 
b OSARGs. but the possibility of self-excitation cannot be ruled 
out. The argument for of the latter interpretation is that peaks are 
narrow, suggesting long-lived modes. 

In the lower panel we show the frequency-luminosity de- 
pendence for the first and second overtones along the three 
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Fig. 2. The b[ and sequences in the v-W/ plane compared with 
the expected positions of power calculated with Eq.(2) {top) and 
with frequencies of the first two radial overtones, V2 and V3 (bot- 
tom). The dots show all significant peaks at the S/N=5 level. 
Some objects yield points on both sequences. 



isochrones. The frequencies were calculated with our nonadi- 
abatic pulsation code (Dziembowski, 1977) for envelope mod- 
els with surface parameters from the isochronesQThe data sug- 
gest that b2 and b3 sequences may be connected with excitation 
of these two overtones. The frequency of the specified mode 
is roughly proportional to the square of mean density, hence 
Vk oc 7T 3 ' 4 . while v max oc L. Therefore, with increasing lumi- 
nosity, the order of modes preferred by SLOs decreases. 

Another characteristic of the SLOs is the dependence of os- 
cillation amplitudes on stellar parameters. Kjeldsen & Bedding 
(1995) derived a semi-empirical expression for the bolometric 
luminosity amplitude of 6L oc (LIM)T~Z. Samadi et al. (2007) 
considered a more general dependence of amplitude on the 
luminosity-to-mass ratio, (L/M) s . They show that the s-exponent 
(sv in their notation) is determined by the eddy time-correlation 
form. The Gaussian form yields s = 1, while the Lorentzian 
s = 0.7. They also show that the latter value fits the ground- 
based data on radial velocity amplitudes better. The space data 
on luminosity amplitudes yield discrepant conclusions. Stello et 
al. (2010) find s « 0.7 from the Kepler data, while Mosser et al. 
(2010) found s = 0.89 + 0.02 from their analysis of the C0R0T 
data. 



1 The parameters do not include the surface helium abundance, which 
in the RGB stars may be up to ten percent higher than the initial one 
because of the dredge-up of the nuclear processed gas. We checked that 
such an increase has no visible effect on the plots shown in this paper. 
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Fig. 3. The highet amplitudes of the peaks on the b2 or b3 se- 
quences as a function of the reddening-free magnitude Wj com- 
pared with the expected values for the solar-like oscillations cal- 
culated with Eq. (3). In the upper panel, the results of calcula- 
tions with the three indicated values of s for models along one 
of the three isochrones (Age=4Gy, Z=0.004) are shown. In the 
lower panel, results for the three isochrones at s = 0.93 are com- 
pared 

The OGLE-III Catalog contains data on mode amplitudes in 
the / band. We compared these data with the prediction from a 
generalized form of Eq.(8) of Kjeldsen & Bedding (1995), which 
at X = 800 nm yields 

ft) 1 *™**-* «> 

The plots in Fig. 3 show that the OGLE data follow the relation 
with s « 0.9. Both, s = 0.7 and 1 are excluded. Our result is, 
thus, similar to that of Mosser et al. (2010). 

3. Modes excited in luminous red giants 

Identification of modes responsible for the peaks in power spec- 
tra such as those shown in Fig. 1 would give us a basis for a real 
seismic probing of individual luminous red giants. The plots in 
the bottom panel of Fig. 2 suggest that the dominant peaks are 
associated with the first and second overtones of radial pulsation. 
If it was true, the frequency distance between the peaks would 
be the large separation, Av. In lower luminosity stars, more than 
one large separation is measured. Their mean value, < Av >, is 
regarded as an important seismic parameter. However, the mean 
distance between sequences b2 and b3 at specified W] cannot 



be identified with < Av > because it does not refer to the same 
object. We need the distance between the dominant peaks mea- 
sured separately in individual power spectra. With this in mind, 
we selected over 100 objects with the most significant (S /N > 7) 
peaks within the two sequences. Plots in Fig. 4 show that, al- 
though the observed Wi(v2) relation is reproduced well by the 
models, the difference between the frequencies, which is shown 
in the lower panel, is by some 20 percent less than the calcu- 
lated difference between the second and the first overtones. This 
difference is quite a robust quantity depending almost only on 
v\ . It is only weakly dependent on the isochrone parameters, as 
may be seen in the figure, as well as on details of modeling the 
envelope and its pulsation. 

Our first suspicion has been that the discrepancy may be 
due to inadequacies of our codes, which were written long ago 
(Paczyhski 1969; Dziembowski 1977). The microphysics has 
been updated, but still the turbulent pressure and the Lagrangian 
perturbation of the convective flux are ignored in these codes. 
Fortunately, Xiong & Deng (2004) have provided a functional 
expression for radial mode frequencies calculated with their 
much more advanced treatment of convection for a sequence of 
red giant envelope models. Repeating calculation with our code 
for models with the same surface parameters, we found only 
a small difference, which was much too small to explain the 
discrepancy shown in Fig. 4. Also acceptable departures from 
the mixing length parameter adopted in the evolutionary mod- 
els only have a minor effect on the dependence of the vj - V2 
difference on vi. 

The problem with interpreting the b2 and b3 sequences in 
terms of the first two overtones could have been noticed earlier 
with a closer look at the Petersen Diagrams shown by Soszynski 
et al. (2004). In such diagrams, period ratios of all significant 
peaks detected in individual objects are plotted against the longer 
periods. In Fig. 1 1 of that paper, we may see a prominent cloud 
of points corresponding to the b2 and b3 sequences near P2/P3 ~ 
0.7. However, most of them are located above this value, whereas 
in the same period range, the models yield values that are always 
below 0.7. The mean value for the cloud of 0.73 is not reached 
until V2 ~ 1/iHz. 

We focused on type b OSARGs because they are, as we be- 
lieve, RGB objects, and we lack needed models of AGB stars. 
Nonetheless, it is still instructive to compare frequencies of type 
a OSARGs with the values calculated along the isochrones in 
the same range of the Wj magnitudes. The essential difference 
between AGB and RGB stars concerns the deep interior, which 
has no effect on radial mode frequencies. The former objects 
have helium-exhausted cores, while the latter are either in the 
phase of core helium burning or still ahead of it. In Fig. 5, we 
compare data on type a OSARGs with results of the the same 
calculations as were used for type b. Points in the upper panel 
seem to follow a single isochrone, which should be somewhat 
younger and/or have higher Z than that for Z = 0.008. This sug- 
gests that at least some of type a OSARGs located below TRGB 
are in fact the first ascent objects before helium ignition. The 
data on frequency differences shown in the lower panel reveal 
the same problem with interpretation in terms of the frequency 
difference between the first two radial overtones. The discrep- 
ancy is only somewhat smaller but there is much less spread in 
the data. So the problem is not less severe. 

We regard the discrepancy revealed in Figs. 4 and 5 as seri- 
ous. Any other assignment of radial mode orders makes the dis- 
crepancy worse. If no modification in stellar models can remove 
the problem with interpretation in terms of the radial modes, 
the option is to consider non-radial modes trapped in the en- 
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Fig. 4. Frequencies of the most significant peaks within the b2 
and b3 sequences in selected objects (dots) compared with fre- 
quencies of the first two overtones in the models along the three 
isochrones. The top panel shows the reddening-free magnitude 
Wi as a function of the lower frequency. The bottom one shows 
the frequency difference. 



velope. Such modes may avoid large radiative losses in the inte- 
rior (Dziembowski et al. 2001, Dupret et al 2009). Though it is 
difficult to explain why such modes could be preferentially ex- 
cited, we did consider this possibility. We limited ourselves to 
dipolar modes because they exhibit relatively small amplitude 
reduction caused by cancelation of opposite sign contributions 
to disk-averaged value. The plots in Fig. 6 show that allowing 
dipolar mode interpretation of the dominant peaks does not solve 
the problerrQ. The vast majority of the data points shown in the 
lower panel falls far from any of the three lines representing cal- 
culated frequency differences. 



4. Discussion and conclusions 

Small-amplitude variability is a common phenomenon in all 
stars along the RGB. Data from the ongoing space missions 
CoRoT and Kepler filled the gap that existed between low- 
luminosity objects, whose oscillations were found in radial ve- 
locity data, and ones located near the tip detected by means of 
massive photometry. The small-amplitude variability observed at 
lower luminosity is always interpreted, as in the Sun, as a mani- 
festation of a stochastic excitation of acoustic modes. 



2 The dipolar mode frequencies were calculated with a new version 
of our code for oscillations in unfitted envelope models that does not 
assume the Cowling approximation. 
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Fig. 5. Same as Fig. 4 but for peaks within the a2 and se- 
quences. 



The first question we addressed in this paper was whether 
oscillations observed in the high-luminosity red giants from the 
OGLE-III catalog (OSARGs) follow the pattern found at lower 
luminosity. We showed that for these stars the oscillation power 
is concentrated in the frequency range of the first two overtones, 
which is consistent with extrapolation of the v max (L) relation 
established for less luminous objects. For the solar-like oscil- 
lations, amplitude increase with luminosity is described by the 
A oc U relation, with s in the [0.7-1.0] range. Samadi et al.(2007) 
show that the value s depends on properties of turbulence, which 
drives the oscillations, and that data for low-luminosity objects 
are best fitted with s x 0.7. De Ridder et al. (2009) obtained a 
similar result for more luminous stars. Our value s « 0.9 for the 
OGLE-III sample of RGB stars is similar to the one found by 
Mosser et al. (2010) for fainter objects. These arguments, in ad- 
dition to the apparent similarity of power spectra shown in Fig. 
1 , support the idea that the variability of stars near the tip of RGB 
is caused by a stochastic excitation of acoustic modes. 

We encountered a problem in our effort to assign radial mode 
frequencies to the two dominant peaks in the OSARGs power 
spectra. The observed frequency dependence on luminosity is 
not far from what is expected for first two radial overtones, but 
the frequency difference is by some 20% less than that calculated 
for the models. The calculated frequency difference was found 
to be only weakly dependent on model parameters. We could 
not find any reasonable modification in models that could lead 
us in the right direction. We checked that the frequencies calcu- 
lated with our code are in good agreement with ones calculated 
for stellar envelope models with the same surface parameters by 
Xiong & Deng (2007). We considered the possibility that one or 
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Fig. 6. Same data as in Fig. 4 compared with frequencies of ra- 
dial (I = 0) and dipolar (( — 1) modes in models along the 
Age=4Gy, Z=0.008 isochrone. The second subscript at v corre- 
sponds to the angular degree, ( 



both peaks are stem from the dipolar mode excitation and found 
that this does not solve the problem either. 

Our failure does not mean the dominant peaks in OSARGs' 
power spectra cannot be associated with low-degree modes of 
stellar oscillations. The problem may still be caused by inade- 
quacy of our models. We do encourage further efforts along this 
line because eigenmode frequencies are indeed the most impor- 
tant seismic observables. 

The alternative, which may seem less attractive, is that the 
peaks are not caused by stellar normal mode excitation. Most 
of narrow peaks in red giants' power spectra from the space 
missions have not been associated with stellar modes either. 
However, even if this is so, the data on red giant variability re- 
main of considerable interest for stellar physics. They certainly 
carry information that deserves to be extracted. In the case of 
OSARGs, the data demanding interpretation are the seven ridges 
in the frequency-magnitude plane. 
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